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Abstract. We examine the effect of point source confusion on cluster detection in Sunyaev-Zel'dovich (SZ) surveys. 
A filter matched to the spatial and spectral characteristics of the SZ signal optimally extracts clusters from the 
astrophysical backgrounds. We calculate the expected confusion (point source and primary cosmic microwave 
background [CMB]) noise through this filter and quantify its effect on the detection threshold for both single and 
multiple frequency surveys. Extrapolating current radio counts, we estimate that confusion from sources below 
~ 100 fj,Jy limits single-frequency surveys to la detection thresholds of Y ~ 3 x 10~ 6 arcmin 2 at 30 GHz and 
Y ~ 10~ 5 arcmin 2 at 15 GHz (for unresolved clusters in a 2 arcmin beam); these numbers are highly uncertain, 
and an extrapolation with flatter counts leads to much lower confusion limits. Bolometer surveys must contend 
with an important population of infrared point sources. We find that a three-band matched filter with 1 arcminute 
resolution (in each band) efficiently reduces confusion, but does not eliminate it: residual point source and CMB 
fluctuations contribute significantly the total filter noise. In this light, we find that a 3-band filter with a low- 
frequency channel (e.g, 90+150+220 GHz) extracts clusters more effectively than one with a high frequency 
channel (e.g, 150+220+300 GHz). 



Key words. 
1. Introduction 

Galaxy cluster surveys based on the Sunyaev- 
Zel'dovich (SZ) effect jSunyaev & Zeldovich 197*51 
|Sunyaev fe Zeldovich 1972| for comprehensive reviews , 
see IBirkinshaw 19991 |Carlstrom, Holder fc Reese 2002| > 
will soon supply large, homogeneous catalogs out 
to redshifts well beyond unity IjBarbosa et al. 1 996 
lEke et al. 1*99*6* IColafrancesco et al. 19971 

IHolder et al. 20001 [Bartlett 20011 IKneissl et al. 20011 
|Benson, Reichardt fc Kamionkowski 2002*) . Eagerly 
awaited, these surveys will probe dark energy and 
its evolution, and shed new light on galaxy for- 
mation HHaiman et al - 20011 [Weller fc Battye 2003| 
[Wang et al. 2 004). The instruments designed for these 
observations are of two types: dedicated interferometer 
arrays surveying at a single frequency (e.g., AMI, AMiBA, 
and SZA) and bolometer systems operating over several 
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millimeter bands (e.g., ACBAR, ACT, APEX, Olimpo, 
Planck, SPT) 1 . Ground-based and balloon-borne in- 
struments are expected to find up to several thousands 
of clusters, while the Planck mission will produce an 
essentially all-sky catalog of several 10 4 clusters by the 
end of the decade. 

An important issue facing these surveys is cluster 
detection in the presence of other astrophysical fore- 
grounds/backgrounds. Except for the very nearby ones, 
clusters will appear as extended sources over arcminute 
scales. Power in diffuse Galactic emission (synchrotron, 
free-free and dust emission) and in the primary cosmic 
microwave background (CMB) anisotropy falls on these 
scales and the clusters can be efficiently extracted using 
an adapted spatial filter ([H aehnelt fc Tegmark 1996[ 
IHerranz et al. 20021 ISchafer et al. 2004jl : fluctuations 



A list of experiment web pages is given at the end of the 
reference section. 
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caused by point sources (radio and infrared galaxies), 
on the other hand, are important on these scales and 
represent a potentially serious source of confusion 
for SZ cluster searches (|Knox, Holder fc Church 2004| 



2. Detecting Clusters: matched SZ filters 

The SZ effect is caused by the hot gas (T ~ 1 - 10 keV) 
contained in galaxy clusters known as the intracluster 



I White k Majumdar 2004 | |Aghanim , Hansen fc Lagac hf200l edium ( ICM ) ; electrons in this S as up-scatter CMB pho- 

-tons and create a unique spectral distortion that is nega- 



The two kinds of SZ survey instruments deal with 
this problem in different ways. Single frequency surveys 
must individually identify and remove point sources with 
high angular resolution (better than 1 arcminute) observa- 
tions, which interferometers achieve by incorporating sev- 
eral long baselines in their antenna array. Operating at 
relatively low frequencies (15 GHz for AMI, 30 GHz for 
the SZA and 90 GHz for AMiBA), these surveys will con- 
tend with the radio galaxy population. Bolometer-based 
instruments will not have the angular resolution needed to 
spatially separate point sources from galaxy clusters; they 
must instead rely on spectral information. In their mil- 
limeter wavelength bands, these instruments will contend 
with the poorly known far-infrared point source popula- 
tion. 



In this paper we quantify the effect of point source 
confusion on cluster detection in SZ surveys. We shall 
only consider the effect of the random field population, 
but we note that point sources are expected to pref- 
erentially reside in the clusters themselves, locally rais- 
ing the effective noise level; we leave examination of 
this effect to a future work. White & Majumdar (2004) 
recently calculated the expected confusion due to IR 
point sources as a function of frequency, while Knox 
et al. (2004) and Aghanim et al. (2004) studied their 
influence on the measurement of SZ cluster parame- 
ters. We extend this work by considering cluster de- 
tection with an optimal filter spatially and spectrally 
(for multi-frequency surveys) matched to the thermal SZ 
signal Paehnelt fc Tegmark 19961 IHerranz et al. 20021 
ISchafer etHfl"" 2004J1 . Using the matched filter, we quantify 
the confusion noise induced by extragalactic point sources 
for both single-frequency radio and multi-frequency 
bolometer SZ surveys. 

We begin by briefly describing our matched filter and 
cluster detection routine, leaving details to Melin et al. 
(2005). In Section |21 we present our point source model, 
based on recent number counts in the radio and far- 
infrared. We then calculate the confusion noise through 
the filter to examine its importance for future SZ surveys 
(Section In the last section, we summarize our main 
results and discuss implications for SZ surveying. 



tive at radio wavelengths and positive in the submillime- 
ter, with a zero-crossing near 220 GHz. The form of this 
distortion is universal (in the non-relativistic limit appli- 
cable to most clusters), while the amplitude is given by 
the Compton y parameter, an integral of the gas pres- 
sure along the line-of-sight. In a SZ survey, clusters will 
appear as sources extended over arcminute scales (apart 
from the very nearby objects, which are already known) 
with brightness profile 



Ai„(x) = y{x)j„ 



(1) 



relative to the mean CMB brightness. Here y{x) is the 
Compton y parameter at position x (a 2D vector on the 
sky) and j v is the SZ spectral function evaluated at the 
observation frequency v. 

A SZ survey will produce maps of the sky at one or 
more frequencies 2 . We model the survey data as a vector 
field M(x) whose components correspond to these maps 
(for a single-frequency survey, the data is a scalar field). 
It is a sum of cluster profiles at positions Xi plus noise and 
foregrounds: M(x) = J2i Vi{ x - x i)jv + N(x), where j v 
is the column vector with components given by j„ evalu- 
ated at the different observation frequencies. The vector 
N(x) includes all non-SZ signals as well as instrumen- 
tal noise; we model it as a stationary random variable 
with zero mean: (N)n = 0, where the average is taken 
over realizations of both the instrumental noise and fore- 
ground fields. We thus assume that the mean intensity of 
the map is zero, i.e., that the zero mode has been taken out 
by the observations 3 . Although the model of a stationary 
random variable applies to the primary CMB anisotropy 
and (perhaps) the noise, one may question its suitabil- 
ity for Galactic foregrounds; it does, all the same, seem 
a reasonable approximation for fluctuations around the 
mean foreground intensity over angular scales pertinent 
to galaxy cluster detection (arcminutes). 



Interferometers actually observe visibilities in the Fourier 
plane, although in this paper we will model interferometer data 
by its image-plane map. This should be a reasonable approx- 
imation if sampling in the Fourier plane is sufficiently good. 
However, it should be emphasized that the question of which 
space is best suited for cluster detection is important and cur- 
rently under study. 

3 In principle this affects the measured cluster Y values, but 
in practice any effect is small for upcoming surveys. 
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Fig. 1. Radial profiles of single frequency and multiple frequency matched niters for a cluster of 9 — 1 arcmin. In the 
left-hand panel we see the spatial weighting used by the single frequency filter to optimally extract the cluster from 
the background (radio point sources + CMB) and noise. The filter is arbitrarily normalized to unity at the origin 
and the beam has a 2 arcmin FWHM. The right-hand panel shows the elements of the 3-band filter with 150, 220 
and 300 GHz channels. The three functions are arbitrarily normalized to the peak value of the 150 GHz filter (black 
curve) , and each channel has a 1 arcmin FWHM beam. In this case we see how the filter uses both spatial and spectral 
information to optimally extract the cluster. 



We now wish to use both spatial and spectral infor- 
mation (when available) to best extract clusters from our 
survey. Consider a cluster of characteristic size 9 (in the 
following, we take this to be the core radius of a /3-profile) 
and central y-value y situated at an arbitrary point Xi 
on the sky. We build a filter i bg(x) (in general a column 
vector in frequency space) that returns an estimate, y Q , of 
y Q when centered on the cluster: 



'Jo = J d 2 x * e *(a: - Xi) ■ M(x) 



(2) 



where superscript t indicates a transpose (with complex 
conjugation when necessary). Suppose that clusters of this 
characteristic size are described by an average angular 
form Tg(x), by which we mean: Tg(x) = (y(x)/y )c, 
where the average is over many clusters of size 9. We match 
the filter to this angular form Tg and to the spectral shape 
of the SZ signal, requiring an unbiased estimate of the cen- 
tral y value: (y ) = y Q , where now the average is over both 
total noise and cluster (of size 9) ensembles. The filter is 
then uniquely determined by further demanding a mini- 
mum variance estimate. 

The result expressed in Fourier space (the flat 
sky approximation is reasonable on cluster angular 



scales) is (|Haehnelt fc Tegmark 1996|IHerranz et al. 2002| 
IMelin fit, al. 2005ji 7 ~ 



* fl (fe) =a 2 e P- 1 (k)-F e {k) 
where 

F e (k) = j v T e {k) 



J d 2 kF e \k)-P- 1 -F e (k) 



-1/2 



(3) 

(4) 
(5) 



with P(k) being the noise power spectrum, a matrix 
in frequency space with components Pij defined by 
(Ni(k)Nf(k')) N = Pij(k)5(k - k'). Note that this last 
expression treats, again, the astrophysical foregrounds as 
stationary fields. The quantity ag gives the the total noise 
variance through the filter. We write the noise power spec- 
trum as a sum P tJ = Pf olse ^j + B t {k)B* {k)P*f y , where 
pnoisc re p resen ts the instrumental noise power, B(k) the 
observational beam and P^ y gives the foreground power 
(non-SZ signal) between channels i and j. As written, we 
assume uncorrelated instrumental noise between observa- 
tion frequencies. 

Our aim is to quantify the effect of point source con- 
fusion on cluster detection using this filter. To this end, 
we ignore diffuse Galactic emission, which is small on 
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cluster scales, and only include primary CMB temper- 
ature anisotropy and point source fluctuations in the 
sky power P^ y ■ For our numerical results, we adopt a 
standard flat concordance CMB power spectrum (O^ = 
1 - n M = 0.7, h = 0.7, e.g., IFreeriman et al. 20011 
|Spergel et al. 20031 ) an d employ a cluster template based 
on the spherically symmetric /3-model with core radius 9 
and p = 2/3: T e (x) = (1 + \x\ 2 / O 2 )-^ 3 -^/ 2 . 

Two examples of the matched filter for 9 = 1 arcmin 
are given in Fig. ^ one for a single frequency survey with 
a 2 arcmin beam (left-hand panel) and the other for a 3- 
band filter with 1 arcmin beams at 150, 220 and 300 GHz 
(right-hand panel). The filters are circularly symmetric 
because we have chosen a spherical cluster model, and the 
figure shows their radial profiles. We clearly see the spa- 
tial weighting used by the single frequency filter to opti- 
mally extract the cluster from the noise and from the point 
source and CMB backgrounds. The multiple frequency fil- 
ter \&e is a 3-element column vector containing filters for 
each individual frequency Their radial profiles are shown 
in the right-hand panel arbitrarily normalized to the peak 
of the 150 GHz filter. The map from each band is filtered 
by the corresponding function and the results are then 
summed to produce the final filter output. We see here 
how the filter uses both spectral and spatial weighting to 
optimally extract the cluster signal. 



3. Point Sources 

Two different extragalactic point source populations af- 
fect SZ observations (see Figure EJ- At frequencies be- 
low ~ 100 GHz, radio galaxies and quasars dominate 
the source counts, while at higher frequencies dusty IR 
galaxies become more important. The spectral depen- 
dence of source flux density in both populations is of- 
ten modeled as a power law, S v cx v a , with spec- 
tral index a. Radio sources tend to have falling spectra 
with a < perbig fc Readhead 1992| |Taylor et al. 2001| 
IMason et al. 2003|) . but flat and inverted spectra with 
a > appear more prominent with observing frequency 
^Bennett et al. 20031 ITTushkin 20 03 ). At millimeter wave- 
lengths we observe the dust emission of IR galaxies in the 
Rayleigh- Jeans with rising spectra characterized by a ~ 
3-4 HVlahakis, Dunne fc Eales 2004"| ) . As will be seen, we 
need information on these extragalactic sources down to 
mjy flux densities and below; unfortunately, neither the 
distribution of a nor the source counts are well known for 
cither population at these flux densities and at frequencies 
of interest for SZ observations (~ 10 — 300 GHz). 
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Fig. 2. Integrated radio and IR source counts. The lower 
solid red line results from the differential counts of Eq. © , 
while the dashed red line corresponds to the model of 
Eq. J7J). Measured counts at 30 GHz from CBI and the 
VSA are shown as the hashed blue boxes. The upper solid 
red line gives the submillimeter source counts from Eq. (JSJ . 
The diamonds show the measured counts at 850 ^m from 
the HDF-North Porys et al. 2003|) and the stars those 
from the 8mJv-survev ljScoT^^tT"2002() . 

3.1. Radio Galaxies 

Bennett et al. (2003) summarize the radio source counts 
at 30 - 40 G Hz from WMAP (5„ > 1 Jy), DASI 
(S„ > 100 mJyIKovac et al. 200211, V SA (S„ > 50 mjy, 
ITaylor et al. "20031 ICleary et al. 2004} and CBI (S u > 
10 mJv. TMlisoT^t£Ll^20'03|l . Fitting to these data, Knox 
et al. (2004) find 



dN 






dS v 


r So 





(6) 



with iV = 30 deg~ 2 , S = 1 mjy and 7 = -2.0. The inte- 
grated counts obtained from Eq. JSJ) are shown in Figure^ 
as the lower solid red line, together with the observations 
from CBI and the VSA. The model fits the CBI counts 
well, but lies high relative to the VSA counts at the bright 
end; as shown in Bennett et al. (2003), the counts in fact 
steepen toward higher flux densities, so a pure power law 
only matches the data over a limited range. 

We will be interested in the counts at flux densities 
near 100 fj,3y for calculating the expected confusion noise 
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in upcoming SZ surveys. This requires an extrapolation 
of the observed counts using Eq. © to much lower flux 
densities, which we view with caution. To get a handle 
on the uncertainty associated with this extrapolation, we 
consider an alternate model with a flattening slope toward 
the faint end: 



dN 
dS u 



So V 



(7) 



with the same values of N a , S Q and 7. These counts are 
shown as the red dashed line in Figure [3 which clearly 
provides an equally satisfactory fit to the observations. 
Comparing the confusion noise in the two models will give 
us a sense of the uncertainty in our estimates. 

We adopt Eq. © , alternatively Eq. (JJJ , for the counts 
at v = 30 GHz. Extrapolation to other frequencies suffers 
from uncertainty in the spectral index a of the emission 
law. Typically negative, determinations of a spread over a 
wide range, including positive values. Mason et al. (2003), 
for example, find an average (a) — —0.45, between 1.4 and 
30 GHz, with a dispersion a a = 0.37, which is roughly 
consistent with the observations of Taylor et al. (2001) 
between 1.4 and 15 GHz. The brighter sources seen over 
the higher frequency WMAP bands, on the other hand, 
show much flatter spectra, with a distribution centered on 
a — and a dispersion of a a ~ 0.3 IjBennett et al. 2 003 
sec also lTrushkin 20030 . There is a pressing need for bet- 
ter understanding of the radio source population at CMB 
frequencies. 

3.2. IR Galaxies 

Dusty infrared luminous galaxies dominate the source 
counts at frequencies near 100 GHz and higher. The dust 
in these galaxies is typically heated to temperatures of sev- 
eral tens of Kelvin by their interstellar radiation field. Its 
emission can be characterized with a blackbody spectrum 
modified by a power-law emissivity: S v oc B v (T), where 
~ 1 — 2. In the Rayleigh- Jeans this leads to a steeply 
rising power-law with spectral index a ~ 3—4, from which 
we see that source confusion from the IR population rises 
rapidly with frequency, the implications of which were re- 
cently discussed by White & Majumdar (2004). 

Blank field counts around 10 mJy at 850 /jm where 
obtained by Scott et al. (2002) and Borys et al. (2003) 
using the SCUBA instrument. As discussed by the latter 
authors, the counts are well described by a double power- 
law: 



dN 
dS v 



No 

So 



with N D = 1.5 x 10 4 deg" 2 and S a = 1.8 mJy. This model 
is shown in Figure|21as the upper solid red line, along with 
data taken from the two surveys. 

The SCUBA Local Univers e Galaxy Survey (SLUGS, 
|Vlahakis, Dunne fc Eales 2004| ) finds a broad distribution 
for the dust emissivity index with (/?) ~ 1 and a dispersion 
we take to be 073 ~ 0.2. According to the SLUGS observa- 
tions, optical galaxies tend to have lower spectral indexes 
than IRAS-selected objects; we eye-balled the above num- 
bers to be representative of the population as a whole. 

4. Point Source Confusion 

Point source confusion is caused by random fluctuations 
in the number of unresolved sources in the filter. We now 
study the contribution of this confusion to the overall filter 
noise, erg, as a function of filter scale, 9. As mentioned 
above, we only consider uncorrelated instrumental noise, 
primary CMB anisotropy and point source terms to the 
power spectrum matrix P, whose off-diagonal elements 
are then just sky terms (to be multiplied by the beam): 



,ps 
ij 



(9) 



where Ci is the CMB temperature angular power spec- 
trum and P^ s is the point source term. We quote power 
in units of CMB temperature equivalent and ignore spa- 
tial correlations of the point sources, which means that 
P? s is independent of k. 

To calculate the point source terms P? s , we adopt the 
counts of Eq Jfjjl, alternatively Eq Q, at v r = 30 GHz 
and those of Eq JSJ at i>„ = 350 GHz (850 /urn). Unless 
otherwise specified, spectral indexes follow Gaussian dis- 
tributions with ((a) = 0, a a = 0.3) for radio sources, and 
((a) = 3, a a = 0.2) for IR galaxies (see previous section). 
Then we have 



PI 



where l|Scheuer 1957MCondon 1974jl 



0,ir) 



V dT 



dS„ sl d dS ^ 



(11) 



(r,ir) 



gives the sky temperature variance due to radio or IR 
sources at the fiducial frequencies v r and u lr ] conversion 
to CMB temperature units is made with the Planck func- 
tion B v at T cmb = 2.725 K l|Mather et al. 19990 . This also 
appears in the R-factors: 



(8) R 



(r , ir) _ {dB v JdT)-HdB U] /dTY 



(dB^/dT)- 



(12) 



(i 
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Fig. 3. Noise (la) on the integrated Compton Y parameter as a function of filter scale. We assume 10 /iK noise/beam 
and adopt 1 arcmin beams in all bands, except for the green 3-dot-dashed curve in the right-hand panel. At these 
frequencies, we only include IR sources below Su m = lOOmJy, following the counts of Eq. JHJ), and with a spectral index 
a = 3 and no dispersion. Left Panel: The black dotted line shows the sensitivity at 150 GHz, where the instrumental 
noise and point source confusion are comparable (see Figure 0J). Confusion is effectively eliminated by a 2-band filter, 
as demonstrated by the red dot -dashed line lying just above the pure instrumental noise limit (black 3 dot-dashed 
curve). The upper solid green line shows the result when adding primary CMB anisotropy, in which case the 2-band 
filter is unable to cleanly separate the three astrophysical components. Addition of a third band at 300 GHz (blue 
dashed line) does not greatly improve the result, because both the 220 and 300 GHz channels are dominated by point 
sources, as shown in Figure 01 (see text). Right Panel: Artificially removing the point sources from the 300 GHz 
channel, we see (black dotted line) that the filter gains enough leverage on the three astrophysical signals to drop its 
noise to the instrumental limit (red dot-dashed line). Source confusion being greatly reduced at lower frequencies, a 
3-band filter with 90 GHz (solid green line) performs better; this remains true even when we degrade the 90 GHz 
beam to a more realistic 2 arcmin FWHM (green 3-dot-dashed curve). 



The function G accounts for spectral variations: 

G {T , ir) {x) = J da N (l , il} (a) x 2a (13) 

with A^( r ir ) being the normal distribution for radio and IR 
spectral indexes, respectively. These results assume that 
there is no correlation between flux density S v and spec- 
tral index, a, so that the joint distribution is separable. 

4.1. Single Frequency Surveys 

Single frequency surveys can only use spatial informa- 
tion to control point source confusion. Operating at v — 



15 GHz and 30 GHz, respectively, AMI and SZA will con- 
tend primarily with the radio source population; AMiBA, 
on the other hand, must deal with both radio and IR 
sources at v = 90 GHz. The former two interferometers in- 
clude long baselines dedicated to identifying and removing 
point sources at high spatial frequency on the sky, where 
they are cleanly separated from the more extended cluster 
SZ emission. Point source confusion is then caused by the 
residual population below the subtraction threshold, Sn m . 

Apart from the beam convolution, point source con- 
fusion contributes to the filter variance <jg in the same 
manner as instrumental noise. Using the counts of Eq. Jfjjl 
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we find for the confusion power: 



30,30 



of = (12 jjK ■ arcmin) 



100 fiJy 



P?li 5 = ^-^15,15 -Gr(l/2) 
> (48 fiK • arcmin) 



2 | ^lim 



V100 ^Jy 



(14) 



(15) 



where the numerical value in the last line assumes 
G r (l/2) = 1, while it most certainly is larger. These val- 
ues are comparable to target sensitivities for interferome- 
ter surveys, implying that these instruments will have to 
subtract sources down to ~ 100 /iJy or better. Assuming 
a 0fwhm ~ 2 arcmin synthesized beam, we calculate the 
confusion limit on the integrated Compton Y parame- 
ter for unresolved clusters. We conclude that confusion 
from sources below 100 /iJy will limit SZ sensitivity to 
Y ~ 3 x 10" 6 arcmin 2 at 30 GHz and Y ~ 10~ 5 arcmin 2 
at 15 GHz (both are at la). The former limit is a fac- 
tor ~ 5 greater than the confusion expected from primary 
CMB anisotropy in the concordence model for unresolved 
clusters | |Melin, Bartlett fc Delabrouille 2005| ). 

Adopting the alternative model of Eq. Q) , these con- 
fusion limits drop by a factor ~ 20, and the dependence 
of the variance on the source subtraction threshold ap- 
proaches cx Sf im . There is clearly a large uncertainty asso- 
ciated with extrapolation of the counts to faint flux den- 
sities. In this light, note that as long as the counts do 
not steepen toward lower flux densities, the confusion es- 
timates are, fortunately, dominated by the counts at S\i m , 
rather than at some unknown cut-off at even fainter levels. 
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Fig. 4. Power spectra of the primary CMB anisotropy 
(solid black line) , point source confusion in different bands 
(as labeled) and instrumental noise; the latter corresponds 
to 10 fiK noise/ 1 arcmin lobe, which we take here to be 
the same in all bands. Only IR sources with S v < Sn m = 
100 mJy are included, according to the counts of Eq. JSJ), 
and with a — 3 and no dispersion. The figure shows the 
that the 300 GHz channel is dominated by sources, ex- 
plaining why the 3-band filter with 90 GHz performs more 
efficiently than the one with 300 GHz (see Figure 0). 



4.2. Multiple Frequency Surveys 

Planned bolometer-based surveys will operate over sev- 
eral millimeter and submillimeter bands, allowing them 
to use spectral information to extract clusters from the 
foregrounds. This will in fact be their only way to reduce 
the effects of point source confusion, because they will not 
have the spatial resolution 4 needed for subtracting point 
sources from cluster images. In the multi-frequency case, 
we refer to the optimal spatio-spectral filter as a multi- 
filter. 

Specifically, the multi-filter performs a weighted sum 
designed to remove foregrounds from the SZ cluster signal, 
as illustrated in the right-hand panel of Figure^ Figurc^l 
helps to understand the filter's workings and it will al- 
low us to draw some important conclusions. The Figure 
shows the filter noise ay (lc) in terms of the integrated 

4 Most planned surveys have arcminute resolution. 



Compton Y parameter as a function of filter scale for vari- 
ous band and foreground combinations. This is calculated 
from the filter variance as ay — ce J d 2 x Tg(x). We take 
as representative of planned observations a survey with 
10 /zK instrumental noise per 1 arcmin lobe (FWHM) in 
all bands. We only include the IR source population (dom- 
inant at these frequencies) below Sn m = 100 mJy, assum- 
ing brighter sources are explicitly subtracted; this is well 
above the knee in the differential counts of Eq. ©. We 
further fix, for the present, the spectral index to a = 3 
with zero dispersion (G = 1). 

Consider first the case with just instrumental noise, 
point sources (no CMB) and two frequencies, one at 
v = 150 GHz and the other at the thermal SZ null, 
v = 220 GHz. Point source confusion is severe in each 
individual band, as illustrated by Figure 0J The dotted 
line in the left panel of Figure [3] shows the total filter 
noise with just the 150 GHz band. When both bands are 
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used in the filter, the filter noise ay drops considerably 
(red dot-dashed line), approaching the pure instrumen- 
tal noise limit for the 2-band filter, shown as the black 
dashed-3 dotted line. It is straightforward to show un- 
der these circumstances that the filter performs a direct 
subtraction of the point source signal from the 150 GHz 
channel by extrapolation of the 220 GHz signal using the 
known spectral index a. What we see here is that the 
subtraction is perfect to the instrumental noise limit. 

The sky of course also includes the CMB signal and 
other foregrounds, and the point source spectral index has 
a non-zero distribution, both of which complicate the sub- 
traction. We discuss the second effect below and now add 
CMB anisotropy, keeping a fixed. The 2-band filter is no 
longer able to separately determine the three sky signals 
(SZ, point sources and CMB), with as a consequence a 
significant rise in total SZ noise, as shown by the solid 
green line. Surprisingly, the situation does not improve 
even when we include more information with a third ob- 
serving band at v — 300 GHz (blue dashed line). 

This interesting result is due to the fact that both the 
220 and 300 GHz bands are dominated by point source 
confusion (see Figure 0} - neither provides good infor- 
mation the CMB anisotropy, so the filter remains unable 
to completely separate the three sky signals. We can test 
this conclusion by artificially removing point sources from 
the new channel at 300 GHz. The result is shown in the 
right-hand panel of the Figure as the black dotted line; 
the total SZ noise try drops to a level comparable to the 
level induced just by instrumental noise (red dot-dashed 
line) , indicating that once again the subtraction is almost 
perfect. 

This has important consequences for SZ surveying. 
Observing at high frequencies, such as 300 GHz, is very 
difficult from the ground due to atmospheric effects; mov- 
ing up in frequency, one approaches strong water vapor 
lines. What we have just seen suggests that including 
bands beyond the thermal SZ null may not be worth the 
effort, at least not for detecting clusters. 

To further explore this issue, we replace the 300 GHz 
channel by a 90 GHz band. At this lower frequency, point 
source confusion is greatly reduced and gives the filter a 
better handle on CMB anisotropy; in fact, as shown in 
Figure 01 point source confusion is well below the instru- 
mental noise level at 90 GHz. The green solid line in the 
right-hand panel of Figure [3] shows the new result: this 
three-band filter with 90 GHz performs significantly bet- 
ter than the one with 300 GHz. It is, however, unlikely 
that all three bands will have the same beam size, which 
we have taken as 1 arcmin throughout this discussion. The 
green 3-dot-dashed curve in the right-hand panel of the 



Figure shows the result for a three-band filter with a 2 ar- 
cmin beam at 90 GHz. Even with the larger beam at the 
lower frequency, the result remains qualitatively the same 
- the filter with 90 GHz performs better than the one with 
300 GHz. 

Multi-frequency surveys will certainly include 150 and 
220 GHz bands to cover the maximum decrement and null 
of the thermal SZ signal. We conclude here that a 90 GHz 
band is a more valuable addition than one at 300 GHz for 
cluster extraction, despite a loss in angular resolution at 
the lower frequency. 

As a final note, we consider the effects of dispersion 
in the source spectral index a. With dispersion, the fil- 
ter can no longer perform a perfect subtraction by ex- 
trapolation across bands; it must instead find appropriate 
frequency weights for a statistically optimal subtraction. 
In the point source power spectrum matrix P^ s , only the 
self power v . remains unaffected; other diagonal ele- 
ments will increase, while correlations between bands (in 
the off-diagonal elements) decrease. We therefore expect 
the filter's performance to decline. 

We examined the importance of this effect using 
the function G\ x defined in Eq. (|13fl . Dispersions up to 
<j a — 0.5 increase matrix elements involving the 150 
and 220 GHz bands of the power spectrum matrix by 
at most 20%, relative to their values with zero disper- 
sion. The 90 GHz channel is of course the most affected: 
the auto-power element -P^ S g increases by more than 
50% for the same dispersion. Nevertheless, when running 
the filter combinations of Figure |3| we find only a small 
change in the filter noise curves, barely perceptible by eye. 
We conclude that even rather large variations in the fre- 
quency dependence of individual source spectra does not 
significantly increase confusion noise through the filter. In 
more general terms, this also suggests that our confusion 
estimates are not strongly dependent on the foreground 
model. 

5. Discussion and Conclusion 

Primary CMB anisotropy and extragalactic point sources 
are the most important foregrounds for SZ surveys. Point 
source confusion is a particularly critical issue because it 
rises rapidly on cluster angular scales. We have quantified 
its importance for both single frequency and multiple fre- 
quency surveys using current estimates of the radio and 
IR source counts and an optimal matched multi-filter for 
cluster extraction. Our main conclusion are: 

— The expected confusion level from radio point sources 
is uncertain due to lack of information on the 
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counts at required flux densities: a power-law ex- 
trapolation of the observed counts (Eq. 0) leads 
to a confusion-limited SZ sensitivity at 30 GHz of 
7~3x 10 -6 arcmin 2 from unresolved sources below 
Sum = 100 /iJy (la for clusters unresolved by a 2 ar- 
cmin beam); at 15 GHz the corresponding sensitivity is 
Y ~ 10 -5 arcmin 2 . An alternative model in which the 
counts flatten at low flux density (Eq.|7J| predicts much 
lower confusion limits, reduced by a factor of ~ 20, in 
which case CMB confusion becomes the limiting factor 
flMelin, Bartlett fc Delabrouille 2005| ). 

— Currently planned multiple frequency bolometer sur- 
veys (e.g., 10 /iK instrumental noise/ 1 arcmin beam) 
greatly reduce the large confusion noise from IR 
sources and can attain 1 a sensitivity of F ~ 7 x 
10~ 6 arcmin 2 for unresolved clusters (see Figure [^J. 
Nevertheless, both residual point source and CMB 
anisotropy confusion significantly increase the total fil- 
ter noise, degrading survey sensitivity from the ulti- 
mate instrumental noise limit. 

— For multiple frequency surveys, observation bands be- 
low the thermal SZ null (e.g., 90+150+220 GHz) are 
more efficient for cluster extraction than observations 
at higher frequencies (e.g., 150+220+300 GHz), de- 
spite a potential loss of angular resolution in the 
90 GHz channel (e.g., 2 arcmin beam instead of our 
fiducial 1 arcmin beam). 

The minimum detectable cluster mass as a function 
of redshift, Mdet(z), is a key characteristic of a survey. 
To further illustrate the last point, we show Mdet(z) in 
Figure |S] for the 3-band multi-filter. In making this fig- 
ure, we adopted simple self-similar relations for Y(M, z) 
and 6> C (M, z) flMelin, Bartlett fc Delabrouille 2005| ); given 
the potentially important theoretical uncertainty in these 
relations, the absolute positioning of the curves should be 
viewed with caution - more robust and more pertinent 
to our discussion are their relative positions. The 3-band 
filter with 90 GHz (black dashed line) gains mass sen- 
sitivity compared to the filter with 300 GHz (solid blue 
line), assuming 1 arcmin beams in all bands. It does not, 
however, reach the ideal noise limit (red dot-dashed line), 
due to residual point source and CMB confusion through 
the filter. A 3-band filter with lower angular resolution 
at 90 GHz (2 arcmin beam), a more probable observing 
situation, does somewhat worse, but still better than the 
3-band filter with 300 GHz (and 1 arcmin beams in all 
bands). This is an important consideration given the dif- 
ficulty imposed by atmospheric effects on observations at 
high frequencies. 

With the S/N > 5 detection threshold [i.e., 5ery(0)], 
we find 15 clusters/deg 2 for the filter with 300 GHz (1 ar- 



CMB+noise+sources 150+220+300GHz 
CMB+noise+sources 90+150+220GHz 
idem but 2 arcmin beam for 90GHz 
noise 90+150+220GH: 




Fig. 5. Minimum detectable mass as a function of red- 
shift. The upper solid blue curve shows the result for the 
3-band filter with a 300 GHz channel, while the middle 
dashed black curve gives the result when this channel is 
replaced by a 90 GHz band. For reference, the lower red 
dash-dotted curve gives the ideal detector-noise limited 
detection mass for the 3-band filter with 90 GHz. The 
3-band filter with 90 GHz shows its higher sensitivity 
(see Figure [3J| , but still suffers from residual foreground 
(CMB) contamination. These results follow for 1 arcmin 
FWHM beams in all bands. For comparison, the black 3- 
dot-dashed line gives the mass limit for the 3-band filter 
with a 90GHz beam of 2 arcmin FWHM, which continues 
to perform better than the 3-band filter with 300 GHz de- 
spite the loss of angular resolution at the lower frequency. 



cmin all bands), 30 clusters/deg 2 for the filter with 90 GHz 
and 1 arcmin beams, and 18 clusters/deg 2 for the 3-band 
filter with a 2 arcmin beam at 90 GHz. There are 47 de- 
tected clusters/deg 2 at the ultimate noise-limit of the 3- 
band filter with a 1 arcmin beam at 90 GHz. These num- 
bers should be compared to the 85 clusters/deg 2 with mass 
> 10 14 solar masses with our model parameters. 

An important consequence of our results is that SZ 
survey selection functions are affected by residual astro- 
physical confusion and are not uniquely determined by 
instrumental properties. Specifically, we have seen that 
even a 3-band bolometer survey with good angular reso- 
lution and optimal filter cluster extraction experiences a 
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mixture of residual point source and primary CMB con- 
fusion. Cluster catalog construction will therefore suffer 
from uncertaintly in astrophysical foreground modeling 
QMelin, Bartlett fc Delabrouille 2005| >. 

In conclusion, our results support the expectation that 
future ground-based SZ surveys will provide rich cluster 
catalogs for cosmological studies. Confusion from point 
sources and primary CMB anisotropy can be greatly re- 
duced by multi-frequency bolometer surveys, but some 
residual point source and CMB anistropy confusion noise 
will affect cluster detection and catalog construction. 
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